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A  simple  model  of  the  interaction  of  the  equatorial  ionosphere  with  the  eastward  F 
region  neutral  wind  in  the  presence  of  evolving  spread  F  bubbles  or  plumes  is  given.  A 
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observed  by  coherent  backscatter  radar  measurements  of  field  aligned  small  scale  irreg¬ 
ularities.  The  essence  of  the  model  is  that  the  plasma  away  firom  the  equatorial  plane 
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ECURITV  CL  ASSiriCATlON  O'  This  r*ct  rwfiwi  Oal*  En/»>4  _  _ 


yj20.  abstract  rCondnuscO 

— (e.g.,  a  background  nighttime  E  region  at  higher  latitudes)  makes  a  finite  contribution  to 
the  magnetic  field  line  integrated  Pedersen  conductivity,  causing  an  incomplete  coupling 
of  the  plasma  motion  to  the  neutral  wind.  The  degree  of  coupling  is  then  a  function  of 
the  Pedersen  conductivities  both  near  the  equator  in  the  F  r^on  and  in  the  hi^er  lat¬ 
itude  E  region,  giving  rise  to  vertical  shears  of  east-west  plasma  motion  having  opposing 
signs  on  either  side  of  the  equatorial  Pedersen  conductivity  Tnayimum.  Evolving  spread  F 
plumes  are  caught  up  in  this  shear  as  they  rise  vertically,  resulting  in  the  characteristic 
“C”  shape  seen  by  backscatter  radar,  and  in  the  westwmd  motion  of  plasma  bubbles 
observed  by  satellite  in  situ  measurements.  Numerical  simulations,  incorporating  an  east¬ 
ward  neutral  wind  in  the  equatorial  F  region  and  E  region  Pedersen  conductivity  effects, 
are  presented  to  further  support  the  model  and  analysis.  The  simulations  also  show  the 
result  that  it  may  be  the  eastward  as  well  as  the  westward  wall  of  the  plume  which  is 
subject  to  secondary  instabilities  in  the  presence  of  an  eastward  neutrd  wind^^  addition, 
even  without  the  neutral  wind,  the  numerical  simulations  show  that  E  regionfPi^ersen 
conductivity  effects  can  result  in  a  slowing  down  of  equatorial  spread  F  and  attendant 
bubble  evolution.  \ 
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NONUNEAR  EQUATORIAL  SPREAF  F: 

THE  EFFECT  OF  NEUTRAL  WINDS  AND 
BACKGROUND  PEDERSEN  CONDUCTIVITY 

1.  Introduction 

In  our  previous  studies  of  evolving  equatorial  spread  F  (ESF)  bubbles 
and  plumes  In  the  equatorial  ionosphere  F  Scannapieco  and  Ossakow,  1976; 
Ossakow  et  al..  1979;  Zalesak  and  Ossakow.  198o3  ,  the  effects  of  the  neutral 
wind  were  neglected.  Rather,  we  focused  our  attention  on  showing  that  the 
motion  and  structure  of  the  experimentally  measured  ESF  environment  (bottom- 
side  and  topside  spread  F,  bubble  formation  and  evolution)  could  be  ex¬ 
plained  in  terms  of  the  nonlinear  evolution  of  the  gravitationally  driven 
collisional  Rayleigh- Taylor  instability.  Through  the  use  of  numerical  simu¬ 
lation  techniques  we  were  able  to  demonstrate  both  qualitative  and  quantita¬ 
tive  agreement  with  the  observations.  We  wished  to  show  that  a  simple  model 
(i.e. ,  using  only  gravity  and  the  bottomslde  background  electron  density 
gradient  as  drivers),  followed  into  the  nonlinear  regime,  could  explain 
observations  that  were  up  to  that  point  inexplicable.  However,  there  are 
some  aspects  of  the  observations  which  we  do  not  see  in  our  previous  simu¬ 
lations. 

First,  there  is  the  tendency  of  ESF  structure  to  drift  eastward  at 
approximately  the  neutral  wind  velocity,  obviously  something  which  could 
not  be  duplicated  in  a  numerical  simulation  which  neglected  neutral  wind 
effects.  Secondly,  there  is  the  curious  tendency  on  the  part  of  radar  back- 
scatter  maps  of  ESF  to  show  plumes  of  backscatter  intensity  which  tilt  east¬ 
ward  with  altitude  at  the  lower  altitudes  and  westward  with  altitude  at  the 
higher  altitudes.  These  structures  were  dubbed  "C's"  and  "fishtails"  by 
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Woodman  and  La  Hoz  [l976]  and  have  also  been  seen  In  the  ALTAIR  backscatter 
maps  of  Tsunoda  [1981] ,  although  Tsunoda  chooses  not  to  regard  the  eastward- 
tilting  and  westward-tilting  structures  as  part  of  the  same  plume.  Addition¬ 
ally,  we  should  point  out  that  McClure  et  al . .  ^19773  have  observed  the 
westward  drift  of  plasma  bubbles  (blte-outs  or  depletions).  We  propose  here 
a  simple  model  of  the  Interaction  of  the  eastward  neutral  wind  at  the  equator 
with  the  equatorial  Ionosphere  which  we  believe  explains  both  of  these  ob¬ 
servations.  At  this  juncture,  we  should  point  out  that  Woodman  and  La  Hoz 
[1976]  ,  Ott  r 1978],  and  Ossakow  and  Chaturvedl  [l978]  hypothesized  that  an 
eastward  neutral  wind  would  produced  westward  drift  of  ESF  bubbles. 

Briefly,  we  find  that  If  the  magnetic  field  line  Integrated  Pedersen 
conductivity  has  a  finite  contribution  from  plasma  which  Is  not  subject  to 
the  equatorial  F  region  neutral  wind  (e.g.,  plasma  at  higher  latitude  E 
regions),  then  the  vertical  polarization  electric  field  driven  by  the  neutral 
wind  at  the  equator  Is  partially  shorted  out  by  this  "background"  E  region 
conductivity,  causing  there  to  be  relative  motion,  or  "slip",  between  the 
plasma  and  the  neutral  wind  at  the  equator.  This  effect  was  first  described 
by  Rlshbeth  Cl97ll  .  Further  Investigation  shows  that  the  degree  of  "slip" 

Is  Inversely  proportional  to  the  "local"  (l.e. ,  equatorial  F  region)  Pedersen 
conductivity,  causing  there  to  be  a  vertical  shear  In  the  plasma  motion  even 
when  there  Is  no  vertical  shear  In  the  neutral  wind.  This  plasma  shear  bends 
any  vertical  structure  about  the  "local"  maximum  In  Pedersen  conductivity, 
giving  rise  to  the  "C's"  and  "fishtails"  seen  on  coherent  backscatter  radar 
maps  (Woodman  and  La  Hoz,  1976;  Tsunoda ,  1981).  Zalesak  et  al.  (1980)  have 
presented  a  preliminary  version  of  this  model. 

In  section  2  we  present  the  geometry  of  the  physical  problem  we  are 
modeling  and  briefly  review  the  relevant  equations  of  motion.  We  also  show 
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that  any  passive  structure  placed  In  the  ambient  equatorial  environment  will, 
in  fact,  be  bent  into  C-shaped  structures.  However,  spread  F  plumes  are 
far  from  passive  structures,  and  it  is  necessary  to  perfora  nonlinear  numeri¬ 
cal  simulations  to  prove  the  case  unequivocally.  These  simulations  are 
presented  in  section  3,  where  we  also  show  the  surprising  result  that,  for 
the  case  studied,  it  may  be  the  eastward  wall  of  the  plume,  as  well  as  the 
westward  wall,  which  is  subject  to  secondary  instabilities.  A  stability 
analysis  which  included  only  the  interaction  of  the  neutral  wind  with  the 
plasma  gradients  in  the  bubble  would  conclude  that  it  should  be  only  the 
westward  wall  of  the  plume  which  is  unstable.  Consideration  of  the  self- 
consistent  polarization  electric  field  of  the  bubble  itself,  as  well  as  of 
gravitational  effects  on  the  tilted  structure  can  cause  the  instability  to 
"switch  sides".  In  section  4  we  present  our  conclusions,  and  in  section  5 
we  discuss  briefly  our  plans  for  future  work. 
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2.  Theory 


In  Figure  1,  we  show  the  geometry  of  the  physical  phenomenon  we  are 
attempting  to  model.  The  equatorial  F  region  plasma  responds  to  the  effects 
of  the  earth's  magnetic  field,  gravity,  collisions  with  the  neutral  atmosphere, 
and  electric  fields.  Since  the  conductivity  along  magnetic  field  lines  Is 
extremely  high,  these  electric  fields  can  depend  on  the  dynamics  of  plasma 
far  from  the  equatorial  region,  but  connected  to  the  equatorial  region  by 
magnetic  field  lines.  We  find  that  the  physical  quantity  dominating  the 
evolution  of  the  colllslonal  Raylelgh-Taylor  Instability  Is  the  magnetic 
field  line  Integrated  Pedersen  conductivity,  and  that  the  primary  contri¬ 
bution  to  that  quantity  comes  from  plasma  In  the  local  region  near  the 
"computational  plane"  shown  In  Figure  1.  This  fact  has  been  the  basis  for 
our  previous  theoretical  and  numerical  studies  of  equatorial  spread  F, 
(Scannapleco  and  Ossakow,  1976;  Ossakow  et  al. .  1979;  Zalesak  and  Ossakow, 

1980)  and  has  enabled  us  to  study  the  phenomena  of  Interest  using  just  a 
single  two-dimensional  computational  plane. 

We  do  not  propose  here  to  analyze  the  problem  In  the  complete  three- 
dimensional  geometry,  but  rather,  as  a  first  step,  to  modify  our  two- 
dimensional  model  to  take  Into  account  the  presence  of  other  plasma,  and 
hence  Pedersen  conductivities  and  forces.  In  regions  far  from  the  equatorial 
plane,  but  connected  to  the  equatorial  F  region  plasma  along  magnetic  field 
lines.  For  Instance  this  could  be  the  northern  and  southern  hemisphere  E 
region  plasma  shown  In  Figure  1.  This  modification  Is  shown  In  Figure  2, 
where  we  show  three  distinct  layers  of  plasma  connected  by  magnetic  field 
lines.  The  center  layer  Is  the  same  computational  plane  as  we  have  used  In 
our  previous  work  (Scannapleco  and  Ossakow.  1976;  Ossakow  et  al..  1979; 
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EQUATORIAL 
SPREAD  F  MODEL 


.  DASHED  LINE  DEPICTS 
E  REGION  PLASMA 


COMPUTATIONAL  PLANE 
'(EQUATORIAL  F  REGION 
PLASMA) 


Fig.  1  —  Diagram  of  the  equatorial  ionosphere  and  of  the  neighboring  regions  which  have 
physical  relevance  to  equatorial  spread  F  (ESF)  processes,  including  the  E  region  plasma 
at  hi^er  and  lower  latitudes.  These  regions  are  electrically  coupled  to  the  equatorial  F 
region  ionosphere  by  the  high  conductivity  along  magnetic  field  lines.  Plasma  is  actually 
distributed  all  along  these  field  lines,  but  in  this  study  we  shaU  make  the  assumption  that 
this  system  can  be  modeled  accurately  by  three  planes  of  plasma  connected  by  straight 
field  lines,  as  shown  in  Figure  2.  One  of  these  thm  layers  (layer  2  in  Figure  2)  is  shown 
here  as  the  “computational  plane.” 
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Fig.  2  —  The  “three  layer”  model  of  the  physical  system  depicted  in  Figure  1.  All  plas¬ 
ma  in  the  vicinity  of  the  equatorial  plane  has  been  compressed  into  layer  2,  while  the 
remaining  northern  and  southern  hemi^here  plasma  has  been  compressed  into  layers  1 
and  3  respectively.  Further,  the  magnetic  field  lines  have  been  strai^tened  so  we  can 
deal  in  cartesian  coordinates  x,  y,  and  z  as  shown  in  the  figure.  The  plasma  in  layers  1 
and  3  is  assumed  to  be  uniform  and  free  of  any  external  driving  force  such  as  a  neutral 
wind.  The  equatorial  layer  2  is  assigned  a  realistic  initial  distribution  of  electron  den¬ 
sity  Ng(y),  and  ion-neutral  collision  frequency,  along  with  a  neutral  wind  which  may 
vary  with  altitude,  but  which  is  taken  to  be  uniform  and  eastward,  and  eqiud  to 
160  m/sec  in  this  study.  Ln  addition,  gravity  points  in  the  negative  y  direction. 
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Zalesak  and  Ossakow,  1980)  and  represents  the  equatorial  nighttime  F  region 
plasma.  The  upper  and  lower  layers  represent  the  remaining  northern  and 


southern  hemisphere  plasma  respectively,  Including  the  E  region  plasma.  The 
problem  Is  still  essentially  two-dimensional  In  that  we  do  not  allow  trans¬ 
port  of  Ions  between  layers,  nor  do  we  allow  any  physical  quantity  to  vary 
with  z  within  a  layer,  where  z  Is  the  direction  along  the  magnetic  field. 

We  do,  however,  allow  electron  currents  to  flow  along  field  lines  between  the 
layers  to  preserve  electrical  neutrality.  Also,  within  the  context  of  this 
model,  we  will  finally  take  the  E  region  layers  to  act  as  a  passive  load, 
l.e.,  we  do  not  allow  for  any  change  In  layers  1  and  3  and  those  layers  are 
assumed  to  remain  uniform.  Thus,  as  a  first  cut  we  are  taking  our  previous 
equatorial  plane  simulations  (Scannapleco  and  Ossakow,  1976;  Ossakow  et  al. , 
1979;  Zalesak  and  Ossakow,  1980)  and  adding  a  passive  E  region  load  to  the 
circuit  to  allow  for  short  circuiting  effects.  Under  the  assumptions  that: 

(1)  the  electric  fields  of  Interest  are  electrostatic  and,  hence,  derivable 
from  a  scalar  potential;  and  (11)  the  conductivity  along  magnetic  field  lines 
Is  extremely  large  and,  hence,  the  potential  Is  constant  along  a  field  line, 
we  are  left  with  a  problem  similar  to  the  multilevel  barium  cloud  striatlon 
problem  f  Lloyd  and  Haerendal ,  1973;  Scannapleco  et  al. ,  1974;  1976;  Doles 
et  al.,  1976].  We  will  now  briefly  derive  the  multilevel  equations,  in 
general  form,  appropriate  to  our  ESF  problem. 

Consider  a  plasma  consisting  of  Ions  and  electrons  Imbedded  In  a  magnetic 
field  aligned  along  the  z  axis.  The  continuity  and  momentum  equations  de¬ 
scribing  the  system  are: 


3n 
_ a 
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where  the  subscript  a  denotes  the  species  (i  for  ions,  e  for  electrons), 
n  is  the  species  number  density,  is  velocity,  v  is  the  recombination  co- 
efficient,  ^  is  the  electric  field,  ^  is  the  gravitational  acceleration,  q 
is  the  species  charge,  is  the  species  collision  frequency  with  the  neutral 
atmosphere,  is  the  neutral  wind  velocity,  c  is  the  speed  of  light,  and 
m  is  the  species  mass.  Note  that  we  have  neglected  finite  temperature  effects 
(pressure),  and  the  effects  of  ion-ion  collisions  and  electron-ion  collisions 
(eventually,  we  will  even  neglect  electron-neutral  collisions).  We  further 
assume  that  we  are  interested  only  in  average  drift  velocities  over  time 
scales  long  compared  to  either  the  mean  time  between  collisions  or  the  gyro- 
period.  In  this  case  we  can  neglect  the  inertial  terms  (the  left  hand  side) 
of  (2),  and  Invert  the  equation  to  obtain  an  algebraic  expression  for  v  : 
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The  vector  subscripts  -i-  and  ||re£er  to  the  components  of  the  vector  which  are 

perpendicular  and  parallel  respectively  to  z.  We  take  "  e  and  ^  -e. 

We  then  assume  that  v  m  0  and  obtain 
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We  now  define  the  perpendicular  current 

i  =  52  n  q  V 
"-x  ^  a  o-ax 
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(16) 

Substituting  (11)  through  (15),  (3)  and  (4)  into  (16), 

and  using  the  quasl- 

neutrality  approximation 
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Eor  our  problem 
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and  we  obtain 
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Since  0.01  £B^<1.0  we  may  neglect  with  respect  to  m^. 

Defining  the  Pedersen  conductivity 

V  . 

_  —  D  /00\ 

ho—  nr 

and  noting  that  ^  ~  ''in^^^l^  obtain 

[m  ro. 

E+-^g+V^—  U 
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Equation  (23)  Isto  be  applied  to  each  of  our  layers  of  plasma.  Referring  to 
Figures  1  and  2  we  see  that  for  layers  1  and  3^  =  0  and  further  that 

o 

j5^  =  -  Dgy  where  g  *  980  cm/sec  and  0<D<1  to  account  for  the  fact  that  £ 

Is  not  perpendicular  to  ^  for  plasma  away  from  the  equatorial  plane  (D  Is 
actually  cos  6^  where  6^  Is  the  dip  angle).  Under  the  assumptions  we  shall 
make  later  It  will  be  seen  that  the  value  of  D  Is  Irrelevant  and  can  be  taken 


to  be  zero.  In  layer  2,  the  equatorial  plane,  we  have  -  g  y,  and  we 

make  the  assumption  that  the  neutral  wind  Is  directed  along  the  x  axis 

(U  >  U  X,  where  x  points  west).  Furthermore,  since  layer  2  Is  taken  to  be 
—  n 

at  F  region  altitudes  where  (Rj^l),  we  can  neglect  In  that  layer 


10 


the  second,  third,  and  fourth  terms  of  equation  (23)  with  respect  to  the  fifth, 
sixth,  and  first,  respectively.  So  we  have  for  the  three  layers: 


where  the  numerical  subscripts  refer  to  the  layer  numbers  depicted  In  Figure  2. 
Quasl-neutrallty  demands  that 
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•  1 
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Integrating  (27)  along  field  lines  and  assuming  that  vanishes  at  z  =  +  “ 
we  obtain 


j  dz  =  0 
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where 
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If  we  model  our  plasma  as  an  array  of  discrete  layers  of  planes  of  plasma  per¬ 
pendicular  to  the  magnetic  field  as  In  Figure  2  we  may  replace  the  Integral  i 

by  a  sum:  | 
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1  ,  Az,  ■  0 

*-Ak  k 


(30) 


where  Az^  is  the  thickness  of  layer  k  measured  along  the  magnetic  field  line. 
By  our  assumption  of  equlpotentlal  magnetic  field  lines  and  electrostatic 
electric  fields 


-xi  “  -J.2  “  -i3 

where  we  have  neglected  the  slight  convergence  of  the  magnetic  field.  Then 
(30)  becomes 


where  the  subscript  b  denotes  the  sum  of  levels  1  and  3  and 
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Also  we  have  defined 
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Note  that  implicit  in  the  above  manipulation  is  the  assumption  that  and 
are  constant  along  a  magnetic  field  line  within  a  given  layer,  as 
we  had  assumed  earlier.  Equations  (1)  and  (32)  constitute  the  system  of 
equations  we  must  solve.  In  general  it  will  be  necessary  to  resort  to 
numerical  means  for  this  task,  but  for  the  case  of  an  unperturbed  laminar 
ionosphere  it  is  both  possible  and  useful  to  find  a  simple  analytic  solution 
to  the  plasma  flow  field,  which  is  an  illuminating  example. 

Suppose  Z  ,  E  ,  S  ,  are  functions  only  of  y  (altitude  in  the 
equatorial  plane).  Then  for  any  set  of  boundary  conditions  on  4*  which  does 
not  itself  impose  an  x-dependence  on  ij>»  we  find  that  4)  =  4>(y)»  Then  (32) 
becomes 


3y 


(X  +  Z  +  2  )  |A=  -  (Z  J2,  —  U  ) 

Pi  P2  P3  3y|  9y  P2  i  e  n' 


(35) 


the  general  solution  of  which  is 
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where  is  a  constant,  and  we  have  dropped  the  subscript  2  on  m^, 


and 
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U^.  Assuming  that  0  as  y  ^  and  demanding  that  (or  equivalently 
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the  total  current)  vanish  at  y  «  +  “  we  get  «  0,  Recalling  that 

-  E  we  obtain 
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The  plasma  motion  produced  by  this  electric  field  Is  given  by 
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f  =  Z  /  (Z  +Z  +Z  ) 
P2  Pi  P2  P3 


(39) 


Note  that  the  plasma  drifts  at  a  fraction  f  of  the  neutral  wind  velocity,  and 

that  that  fraction  Is  simply  the  ratio  of  the  "local"  (l.e.,  equatorial  plane) 

Pedersen  conductivity  to  the  total  field  line  conductivity  on  a  given  field 

line  (Note:  what  we  have  In  mind  here  and  In  the  numerical  simulations  Is 

that  our  magnetic  field  line  Integration  for  the  equatorial  F  region  is  over 

a  few  degrees  In  latitude,  and  that  regions  1  and  3  constitute  the  rest  of 

the  field  line  connected  Ionosphere  as  a  load  on  the  circuit).  This  simple 

equation  has  some  remarkable  consequences  in  terms  of  the  motion  of  structures 

(spread  F  plumes,  for  example)  imbedded  In  the  equatorial  Ionosphere.  Sup- 

pose  that  Z  Is  a  function  of  altitude  with  a  peak  Z  at  altitude  h 
P2  P2 

Suppose  further  that  Z  and  Z  are  constants  such  that  Z  +  Z  *  0.1  z'"*®*. 

Pi  P3  Pi  P3  P2 

and  that  we  Impose  a  uniform  eastward  neutral  wind  of  100  m/sec  on  level  2 
(the  equatorial  plane).  We  now  create  a  model  ionosphere  (see  Table  1)  and 


14 


tabulate  the  eastward  plasma  velocity  as  a  function  of  altitude: 


TABLE  1 


Altitude  (km)  2 

_  -22 _ E2_ 


Eastward  Plasma 
Velocity  (m/sec) 


600 

0.1 

50 

500 

0.5 

83 

400  (h  ) 
max 

1.0 

91 

300 

0.1 

50 

200 

0.01 

9 

Note  that  even  though  there  Is  no  vertical  shear  In  the  neutral  wind  velocity, 

the  plasma  flow  field  contains  a  large  shear  with  opposing  signs  on  either 

side  of  h  .  The  effect  of  this  shear  Is  to  bend  any  passive  vertical 
max 

structure  Imbedded  in  this  flow  field  into  a  "C"  shape  as  depicted  in 

Figure  3  (Also  note  that  for  2  =»  2  =0,  i.e,,  no  E  region,  from  (39) 

PI  P3 

f  «  1  and  the  plasma  moves  at  the  wind  speed  (Rishbeth,  1971)), 

The  above  result  Is  quite  satisfying  In  that  it  offers  a  qualitative 
explanation  of  the  "C's",  "fishtails",  and  other  tilted  structures  seen  by 
Woodman  and  La  Hoz  [l976]  and  Tsunoda  [l98l]  in  their  observations  of  co¬ 
herent  radar  backscatter  from  the  meter-scale  irregularities  associated  with 
ESF  plumes.  However,  the  above  analysis  is  valid  only  for  passive  structures 
imbedded  in  a  laminar  unperturbed  ionosphere,  conditions  which  are  simply 
not  met  in  the  ESF  environment.  Numerical  simulations  are  necessary  to 
prove  the  case  unequivocally. 
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Fig.  3  —  Schematic  diagram  depicting  the'bending  of  a  passive  vertically-aligned 
structure  caught  up  in  a  velocity  shear  pattern  of  the  type  we  believe  exists  in  the 
nighttime  equatorial  ionosphere.  The  neutral  wind  is  eastward  and  uniform  in 
altitude,  and  the  response  of  the  plasma  (depicted  by  arrows)  is  to  move  at  some 
fraction  of  the  neutral  wind  velocity,  that  fraction  being  largest  at  the  altitude 
hmax  maximum  equatorial  plane  Pedersen  conductivity.  The  eastward  plasma 
velocity  falls  off  both  above  and  below  h_„,  as  shown. 


3.  Numerical  Slmulatious 

We  mentioned  before  that  equations  (1)  and  (32)  constitute  the  system 
we  wish  to  solve  numerically.  Let  us  now  be  more  specific.  Equation  (1) 

Is  actually  six  equations  (one  electron  and  one  Ion  equation  for  each  of  our 
three  layers).  By  quasl-neutrallty  (17)  we  can  eliminate  three  of  these  and 
Integrate  either  the  Ion  or  electron  equation  at  each  level;  but  since  we 
have  made  the  assumption  that  v^jj  =  0  (currents  along  field  lines  are  carried 
by  electrons)  we  can  more  conveniently  solve  the  two-dimensional  Ion  continuity 
equation  at  each  layer: 

jit  +  .  („  „ j 

-'k;  k  =  1,2,3  (40) 

In  the  simulations  we  present  here  we  have  set  In  (40)  to  zero  for  simpli¬ 
city  and  because  at  the  F  region  altitudes  we  shall  be  dealing  with,  recombina¬ 
tion  effects  are  negligible. 

We  now  make  one  last  simplifying  assumption:  the  background  E  region 
plasma  (layers  1  and  3)  Is  Initially  uniform  In  density  and  Pedersen  conduc¬ 
tivity,  and  remains  so  during  the  course  of  our  simulation.  This  Is  tantamount 
to  neglecting  compressibility  (Pedersen  mobility)  effects  In  the  E  region 
plasma.  Thus,  we  are  utilizing  layers  1  and  3  as  a  passive  load  In  an  Iono¬ 
spheric  circuit.  In  order  to  allow  for  short  circuiting  effects.  A  true  multi¬ 
level  numerical  code  which  will  model  these  effects  self-conslstently  Is  under 
development.  This  assumption  does  have  hhe  advantage,  though,  of  reducing  (40) 
to  a  single  equation  (since  3n/dt  -  0  for  levels  1  and  3)  and  of  eliminating 
H  and  the  last  two  terms  of  (32)  (since  all  the  terms  subscripted  by  b 
are  constant).  Our  final  pair  of  equations  to  be  solved  numerically  Is  then 


(41) 


V  • 

X 


0 


(42) 


where  all  quantities  except  I  and  I  refer  to  layer  2. 

Pi  P3 

Equation  (41)  Is  solved  numerically  using  the  fully  multl-dimenslonal 
flux-corrected  transport  (FCT)  techniques  of  Zalesak  [l979l.  Briefly, 

FCT  Is  a  technique  originally  developed  by  Boris  and  Book  [l973]  for  solving 
equations  of  the  form  (41)  where  steep  gradients  In  n  are  expected  to  form. 
The  fluxes  used  In  the  algorithm  are  nonlinear  weighted  averages  of  fluxes 
computed  by  high  and  low  order  finite  differences.  The  high  order  fluxes  are 
weighted  as  heavily  as  possible  subject  to  the  constraint  that  nonphysical 
oscillations  are  not  Introduced.  Equation  (42)  Is  solved  using  the  fully 
vectorized  Incomplete  Cholesky  conjugate  gradient  (ICCG)  algorithm  of  Haln 
[l980j,  which  Is  an  extension  of  the  work  of  Kershaw  ^1978^.  This  algorithm 
Is  extremely  fast  and  efficient  for  the  cases  described  below  for  which  the 
neutral  wind  was  set  to  zero;  however,  when  a  finite  eastward  neutral  wind 
was  used  the  ICCG  convergence  rate  became  painfully  slow  and  It  was  necessary 
to  resort  to  the  direct  elliptic  solver  of  Madala  [l978j. 

The  numerical  calculations  to  be  presented  were  performed  on  a  two- 
dimensional  cartesian  mesh  using  40  points  In  the  x  (east-west)  direction  and 
140  points  In  the  y  (vertical  direction) .  The  (uniform)  grid  spacing  was 
3  km  In  the  y  direction,  and  5  km  is  the  x  direction  for  all  calculations. 
Note  that  in  our  previous  work  we  used  2  km  spacing  in  the  y  direction. 

The  bottom  of  the  grid  corresponds  to  253  km  altitude  and  the  top  of  the 
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grid  to  676  km  altitude.  Periodic  boundary  conditions  were  imposed  on  both 
n  and  <j)  in  the  x  direction.  In  the  y  direction  transmissive  boundary  con¬ 
ditions  were  imposed  on  n  (3n/3x  =0),  and  the  normal  derivative  of  (Ji  at  the 
top  boundary  was  chosen  such  that  the  normal  component  of  the  total  current 
(the  sum  over  all  three  layers)  was  zero  there  for  the  unperturbed  state. 

For  the  calculations  with  no  neutral  wind  this  implies  3(J>/3y  =  0  at  the  upper 
boundary.  For  calculations  with  a  neutral  wind  this  implies 


(43) 


at  the  upper  boundary,  where  is  the  Pedersen  conductivity  of  the  initial 
unperturbed  state.  At  the  lower  boundary  we  set  (Ji  =  0  for  all  cases. 

Three  kinds  of  plots  will  be  presented:  (1)  contours  of  constant 
n(x,y,t);  (2)  contours  of  constant  n(x,y,t)/n^(y) ;  and  (3)  contours  of  con¬ 
stant  electrostatic  potential  (()(x,y,t).  Here  n^(y)  is  the  initial  unperturbed 
electron  density  profile  in  layer  2.  Superimposed  on  each  contour  plot  is  a 
dashed  line  depicting  n^(y)  for  reference  purposes.  Our  o^(y)  profile  is 
such  that  the  peak  is  located  at  434  km  altitude,  and  the  minimum  electron 

density  scale  length  L  =  n  (3n  /3y)  ^  is  10  km.  The  ion-neutral  collision 

o  o 

frequency  used  in  the  calculations  is  shown  in  Figure  4.  The  initial 

perturbation  used  to  start  each  calculation  was  a  mode  1  sine  wave  in  the 
X  direction: 


o^*y)  ■  1  -  e  ^  cos  (nx/100) 


(44) 


Three  calculations  were  performed  to  determine  the  effect  of  the 
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ALTITUDE  (km) 


0.01  0.1  I  10  100  200 


ION-NEUTRAL  COLLISION  FREQUENCY.  i;jj^(sec'') 

Fig.  4  —  The  ion-neutral  collision  frequency  (solid  line)  as  a  function 
of  altitude.  The  recombination  coefficient  i'll  was  set  to  zero  for  this 
study  (see  text). 
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background  E  region  plasma  and  of  the  eastward  neutral  wind  on  evolving 
spread  F  bubbles: 

1)  Calculation  2L,  in  which  E  =  E  =  U  =  0.  This  calculation  is 

Pi  P3  n 

identical  to  calculation  2L  of  Zalesak  and  Ossakow  198ol  except 
for  the  difference  in  vertical  grid  spacing  noted  previously. 

2)  Calculation  2LE,  identical  to  calculation  2L  above,  except  that  a 
constant  background  Pedersen  conductivity  has  been  included  such 
that  E  +  Z  =  0.12  E°  ,  where  E°  is  the  maximum  Pedersen 

Pi  P3 

conductivity  in  the  initial  unperturbed  equatorial  plane  (layer  2). 

We  believe  the  value  of  0.12  for  the  relative  background  Pedersen 
conductivity  level  to  be  a  conservative  figure.  Rishbeth  (]l97l] 
used  a  value  of  0.2  as  being  representative  of  nighttime  conditions. 

3)  Calculation  2LEW,  identical  to  calculation  2LE  above,  except  that 
a  uniform  eastward  neutral  wind  of  150  m/sec  was  imposed  over  the 
entire  equatorial  plane  (U^  =  -  150  m/sec).  The  designations 

E  and  W  above  obviously  refer  to  the  presence  of  E  region  plasma  and 
neutral  winds,  respectively. 

Figure  5  shows  isodensity  contours  of  n(x,y,0)  for  our  initial  conditions 
(laminar  ionosphere  ^^(y)  plus  perturbation  (44)).  The  contours  are  labeled 
for  later  reference  purposes.  Note  that  in  this  and  all  subsequent  plots 
we  have  plotted  two  periods  (recall  that  we  have  periodic  boundary  conditions 
in  the  x  direction)  of  the  various  functions.  That  is,  the  40  by  140  mesh 
was  extended  to  80  by  140  for  plotting  purposes  only,  to  facilitate  comparison 
with  plots  of  calculations  run  with  a  neutral  wind,  in  which  structures  move 
across  the  grid. 

Figure  6  shows  Isodensity  contours  of  n(x,y,t)  for  calculation  2L  at 
four  different  times  during  the  simulation.  Figure  7  shows  a  similar  sequence 
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ALTITUDE  (KM) 


EAST  (KM)  VEST 

Fig.  5  —  Iso-electron  density  contours  for  the  initial  perturbed  state  in  layer  2. 

This  represents  the  initial  conditions  for  our  numerical  simulation,  The  con¬ 
tours  are  labeled  in  units  of  electrons/cm^  in  E  format  notation  (l.OEl  =  1  X 
IQl ,  etc.).  The  unperturbed  ionosphere  was  initially  laminar  (independent  of  x, 
the  east-west  direction)  and  is  exhibited  by  the  dashed  curve  showing  N^Cy),  at 
any  point  in  the  east-west  (x)  direction.  This  curve  is  labeled  at  the  top  of  the 
figure.  The  perturbation  has  a  maximum  amplitude  of  e~®  in  relative  electron 
density,  is  a  pure  mode  1  sine  wave  in  x,  and  is  independent  of  altitude  y,  as 
described  in  the  text.  The  observer  is  looking  southward  so  that  B  is  out  of  the 
figure. 


22 


(MY) 


r 


J  1 _ r _ I _ 1 _ I 

-200  -too  0  100  200 

EAST  IKNI  VEST 


i  _ 1 _  •  _ I _ ZJ 

-200  -100  0  100  200 

EAST  IKNI  VEST 


2215  SEC  No  X  I0->CN  * 


EAST  IKNI  VEST 


2560  SEC  No  X  lO'TCN-i 


EAST  IKNI  VEST 


Fig.  7  —  Same  as  Figure  6,  but  for  calculation  2LE  (background  E  region, 
no  wind)  at  1500, 1900,  2215,  and  2560  seconds 


for  calculation  2LE.  (The  reader  may  note  that  the  2L  calculation  of  this 
paper  evolves  at  a  faster  rate  at  late  times  than  the  2L  calculation  of 
Zalesak  and  Ossakow  ^198o3.  This  Is  primarily  due  to  two  Improvements  In  the 
numerical  treatment  of  equations  (41)  and  (42),  Implemented  since  our  previous 
studies:  1)  the  differencing  of  the  Hermltlan  form  (42)  of  the  potential 
equation,  rather  than  the  non-Hermltlan  expansion  we  were  constrained  to  use 
previously,  as  discussed*  In  the  appendix  of  Zalesak  and  Ossakow  fl9803;  and 
2)  Improved  treatment  of  the  continuity  equation  (41)  which  has  enabled  us  to 
further  reduce  the  numerical  diffusion  that  Inevitably  occurs  across  electron 
density  gradients  as  steep  as  those  formed  at  the  edges  of  ESF  plumes  at  late 
times.  We  would  emphas''ze  that  the  conclusions  of  Zalesak  and  Ossakow  [^1980^ 
do  not  depend  on  late-tlme  rise  velocities  and  are,  therefore,  unaffected  by 
this  result).  In  comparing  Figure  6  with  Figure  7  we  are  looking  at  the  ef¬ 
fect  of  a  background  conductivity.  The  most  striking  difference  is  that  of 
the  time  scales,  whereby  2LE  takes  about  70%  more  time  to  achieve  a  600  km 
altitude  plume  than  does  2L.  Qualitatively  this  can  be  understood  In  terms 
of  the  shorting  effect  of  the  background  E  region,  by  which  a  given  current 
can  be  driven  by  a  smaller  electric  field,  which  In  turn  means  smaller  plasma 
velocities.  Almost  as  striking  Is  the  fact  that  the  2LE  plume  bifurcated 
while  the  2L  plume  did  not.  The  inevitability  of  the  bifurcation  in  calcula¬ 
tion  2LE  can  be  seen  even  In  the  very  early  time  plot  at  1500  sec,  where  the 
characteristic  flattening  of  a  significant  number  of  contours  In  the  upper 
portion  of  the  plume,  the  sure  signal  of  Imminent  bifurcation  In  barium  cloud 
studies  r  Zabusky  et  al.*.  1973;  Scannapleco  et  al . ,  1974;  Ossakow  et  al ■ ,  1977; 
McDonald  et  al . ,  198o] ,  can  be  seen.  The  close  similarity  of  the  physics  of 
the  ESF  gravitational  Instability  and  that  of  the  ExB  gradient  drift  In¬ 
stability  associated  with  the  bifurcation  and  strlatlon  process  In  plasma 
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clouds,  has  been  noted  by  Scannapleco  and  Ossakow  [|l976l.  We  shall  draw 
heavily  on  our  knowledge  of  bifurcation  tendencies  in  plasma  clouds 
[* Ossakow  et  al.,  1977;  McDonald  et  al . ,  1980]  when  we  address  the  question  of 
why  the  plume  In  calculation  2LE  bifurcated  while  that  In  calculation  2L  did 
not,  later  In  this  paper. 

For  the  moment  we  note  that  there  are  two  primary  effects  of  the  pre¬ 
sence  of  a  background  conducting  region  (E  region):  1)  electric  fields  every¬ 
where  are  reduced  by  the  shorting  effects  of  the  background  conductivity,  re¬ 
sulting  In  an  overall  slower  evolution  for  the  Instability;  and  2)  electric 
fields  are  reduced  the  most  In  the  regions  where  the  ratio  of  the  equatorial 
plane  conductivity  to  that  of  the  background  plane  Is  smallest,  l.e.,  at  low 
altitudes,  rendering  the  2LE  configuration  Incapable  of  drawing  plasma  from 
extremely  low  altitudes  to  produce  large  depletion  levels  Inside  the  bubble. 
This  can  be  seen  easily  in  comparing  Figures  6  and  7  wherein  we  note  that 
the  isodensity  contours  at  low  altitudes  are  virtually  stationary  in  the  2LE 
case,  whereas  the  2L  configuration  results  In  significant  upward  movement  In 
even  the  lowest  density  plasma  near  the  bottom  of  the  plot.  The  more  effec¬ 
tive  shorting  of  the  electric  fields  by  the  background  layer  In  the  2LE 
case  is  seen  dramatically  In  Fig.  9a  and  9b,  where  we  plot  contours  of 
constant  electrostatic  potential  (p  for  calculations  2L  and  2LE  respectively 
at  early  time.  (The  contour  level  increment  of  (J>  In  this  paper  is  chosen 
such  as  to  divide  the  maximum  excursion  of  from  zero  into  7  equal  Intervals 
The  contours  corresponding  to  positive  values  of  <|)  are  plotted  as  solid  lines 
while  those  corresponding  to  negative  values  are  plotted  as  dashed  lines. 

The  zero  contour  level  Is  suppressed.  For  the  2L  and  2LE  cases  the 
symmetry  of  the  potential  would  cause  the  zero  contour  to  be  simply  two 
vertical  lines.  Since  the  electric  field,  and  hence  plasma  velocity,  is 
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inversely  proportional  to  the  contour  spacing,  this  normalization  allows  us 
to  easily  determine  by  eye  the  rate  at  which  the  upward  velocity  of  plasma 
is  decreasing  with  decreasing  altitude.  It  also  allows  us  to  visualize  the 
global  plasma  flow  field,  since  contours  of  4*  are  essentially  streamlines  of 
the  plasma  flow).  Comparing  Figures  9a  and  9b,  we  note  a  much  more  rapid  de¬ 
crease  in  the  horizontal  component  of  the  electric  field  with  decreasing  al¬ 
titude  in  calculation  2LE  than  in  2L.  The  flow  field  in  2LE  is  mixing  plasma 
over  a  fairly  narrow  altitude  range,  while  that  in  2L  is  drawing  plasma  from 
deep  in  the  ionosphere,  where  the  plasma  densities  are  lowest.  Hence,  we 
should  expect  the  late  time  plume  in  calculation  2L  to  consist  of  plasma  of 
lower  density  (l.e.,  to  have  much  higher  depletion  levels)  than  that  in  cal¬ 
culation  2LE.  That  this  is  Indeed  the  case  can  be  seen  in  Figures  10a  and 
10b,  where  we  compare  Isodenslty  contours  of  n(x,y)/n^(y)  at  late  times  for 
the  two  calculations.  (Contours  of  n/n^  in  this  paper  are  spaced  logarithmi¬ 
cally,  with  solid  lines  representing  depletions  (n/n^<l)  and  dashed  lines  re¬ 
presenting  enhancements  (n/n^>l}.  The  kth  depletion  contour  represents  an 
n/n^  value  of  0.5  ,  while  the  kth  enhancement  contour  represents  an  n/n^ 
value  of  2.0  ).  Although  the  bifurcation  of  the  2LE  plume  makes  the  compari¬ 
son  less  clean  than  it  would  otherwise  be,  it  is  obvious  that  the  depletion 
levels  of  the  2L  plume  are  much  higher  than  those  of  the  2LE  plume.  Depletion 
levels  (1-n/n^)  in  the  upper  portions  of  the  2L  plume  are  greater  than  99.2%, 
while  those  in  the  2LE  plume  are  only  about  94%.  We  conclude  that  the  pre¬ 
sence  of  a  background  conducting  region  results  in  plumes  which  are  both 
slower  to  evolve  and  less  depleted  than  their  no-background  conterparts. 

In  Figure  8  we  present  Isodensity  contours  of  calculation  2LEW  at  times 
similar  to  those  presented  for  the  2LE  calculation.  For  the 
profiles  chosen  we  find  a  peak  in  at  394  km,  40  km  below  the  F2  peak,  with 
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Fig.  9  —  Early  time  contours  of  constant  electrostatic  potential  0(x,y)  for  (a)  calculation 
2L  at  875  sec,  (b)  calculation  2  L£  at  1500  sec  and  (c)  calculation  2  LEW  at  1500  sec. 
The  potential  0o(y)  associated  with  the  unperturbed  initial  conditions  has  been  removed 
in  (c).  The  contour  levd  increment  is  chosen  such  as  to  divide  the  maximum  excursion 
of  0  from  zero  into  7  equal  intervals.  The  contours  corresponding  to  positive  values  of  0 
are  plotted  as  solid  lines,  while  those  corresponding  to  negative  values  are  plotted  as  dashed 
lines.  The  zero  contour  level  is  suppressed. 


falling  off  by  a  factor  of  ten  42  km  belov  and  132  km  above  this  altitude 
P2 

(at  352  km  and  526  km  altitude  respectively).  Using  Equation  (38)  to  approxi¬ 
mate  our  Initial  shear  field  and  using  E  +  Z  =  .12  we  find  eastward 

Pi  P3  2 

plasma  drifts  of  134  m/sec  at  394  km  altitude  and  68  m/sec  at  both  352  and 
526  km  altitude.  The  plasma  shear  is  weaker,  but  over  a  larger  altitude  range, 
above  the  peak  In  Z^  than  below  It.  If  vertical  plasma  plumes  behave  as 
passive  structures,  we  would  expect  a  bending  of  the  structure  around  an 
altitude  of  394  km,  with  a  larger  slope  below  this  altitude  than  above  it. 
Looking  at  Figure  8  we  see  that  this  behavior  is  qualitatively  reproduced,  in 
spite  of  the  fact  that  the  self -cons is teat  polarization  fields  produced  by 
the  plumes  represent  very  large  perturbations  on  the  equilibrium  fields  pro¬ 
ducing  the  plasma  shear.  In  Figure  8  we  have  placed  ourselves  in  a  frame 
moving  at  68  m/sec  eastward  to  minimize  both  computational  errors  and  computer 
time.  In  Figure  10c  we  show  late  time  isodensity  contours  of  n(x,y)/n^(y) 
for  calculation  2LEW,  for  coraparlslon  to  Figures  IQa  and  10b.  The  bending 
of  the  plume  into  a  "C"  shape  about  an  altitude  of  360  km  is  quite  pro¬ 
nounced.  The  fact  that  this  "bending  point"  is  more  than  30  km  below  the 
initial  maximum  in  equatorial  plane  Pedersen  conductivity  is  an  indication  of 
a  nonlinear  interaction  between  plume  rise  and  ambient  plasma  shear.  In  fact, 
this  shift  downward  can  be  understood  qualitatively  as  follows.  The  movement 
of  low  density  plasma  upward  inside  the  plume  is  accompanied  by  the  movement 
of  high  density  plasma  downward  in  the  regions  between  the  plumes  (see  Figures 
9  and  10) .  Since  the  scale  length  over  which  is  decreasing  with  altitude 
is  long  ('^60  km)  compared  to  the  scale  length  over  which  the  electron  density 
is  Increasing  with  altitude  ('v-lO  km)  below  the  F2  peak,  the  effect  of  this 
downward  movement  of  high  density  plasma  is  to  move  the  point  of  maximum 
Pedersen  conductivity  in  the  equatorial  plane  and  hence  the  bending 
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Fig.  10  —  Late  time  contours  of  constant  n{x,y)/ng{y)  for  (a)  calculation  2L  at  1276  sec, 
(b)  calculation  2LE  at  2560  sec,  and  (c)  calculation  2LE  at  2331  sec.  The  contours  are 
spaced  logarithmically,  with  solid  lines  representing  depletions  (n/nQ<  1)  and  dashed  lines 
representing  enhancements  (n/ng>l).  The  depletion  contour  represents  an  n/n^  value 
of  O.S*',  while  the  enhancement  contour  represents  an  n/n^  value  of  2.0*‘ . 


point,  >daw3Afard. 

In  Figure  9  c,  we  show  contours  of  constant  4>  (x,y)-<|>^(y)  at  1500  sec 

for  calculation  2LEW,  for  comparison  to  Figures  9a  and  9b.  Here  ‘P^Cy)  is 

the  initial  equilibrium  electrostatic  potential  of  the  unperturbed  initial 

conditions  (for  the  2L  and  2LE  cases  <fi^  =  0) .  Subtracting  this  quantity 

from  i)>  before  plotting  enables  us  to  examine  the  "underlying"  plume  motion 

on  which  the  shear  associated  with  the  initial  conditions  is  superposed. 

Remarkably,  the  motion  of  plasma  in  the  plume  is  not  purely  upward,  but 

rather  upward  and  westward ,  despite  the  fact  that  we  have  removed  the 

asymmetry-inducing  profile  of  the  equilibrium  initial  wind  field.  Although 

this  simple  analysis  is  crude  (in  that  the  dependence  of  (|>  on  the  plasma 

structure  is  not  linear,  l.e.,  in  Eq.  (42)  if  <|>^(x,y)  and  i<i^(x,y)  are 

solutions  for  (x,y)  and  E®  (x,y)  rerpectlvely,  (t>^(x,y)  +  4'^(x,y)  is  not 

P2  P2 

A  B  • 

a  solution  for  E  (x,y)  +  E  (x,y)),  it  would  seem  to  lend  support  to 
P2  P2 

the  ideas  advanced  by  Woodman  and  La  Hoz  [l9763,  Ossakow  and  Chaturvedi  [1978], 
and  Ott  [l978]  who  proposed  that  a  neutral  wind  whose  eastward  velocity  ex¬ 
ceeded  that  of  the  plasma  would  combine  with  gravity  to  form  an  effective 
gravity  which  pointed  downward  and  eastward,  causing  bubbles  or  plumes  to 
drift  upward  and  westward  relative  to  the  surrounding  plasma.  Thus,  the 
westward  tilt  of  plumes  at  high  altitudes  would  appear  to  be  due  to  both  this 
effect  (since  we  have  shown  that  the  plasma  velocity  does  lag  the  neutral 
wind  velocity)  and  that  of  the  plasma  shear  which  we  have  addressed  earlier. 

We  would  point  out,  however,  that  the  mechanism  of  Woodman  and  La  Hoz  [l976], 
Ossakow  and  Chaturvedi  [l978],  and  Ott  [l978]  cannot  explain  the  eastward 
tilt  of  plumes  with  altitude  at  low  altitudes. 

In  comparing  calculations  2LE  and  2LEW  (Figures  7  and  8),  we  are 
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evaluating  the  effect  of  the  neutral  wind  Itself,  since  the  same  background 
Pedersen  conductivity  was  used  In  both  cases.  We  wish  to  note  the  following 
points:  1)  the  two  calculations  evolve  at  approximately  the  same  rate  In 
time;  2)  the  primary  effect  of  the  wind  Is  to  bend  the  plume  In  calculation 
2LEW  Into  a  "C"  shape,  with  the  upper  part  of  the  "C"  being  much  larger  In 
altitude  extent  and  tilted  markedly  westward;  3)  the  plume  depletion  levels 
are  approximately  the  same  in  both  calculations;  and  4)  the  2LE  plume  bifur¬ 
cated  while  the  2LEW  plume  did  not.  We  shall  address  this  last  question, 
along  with  the  question  of  why  the  plvime  in  calculation  2L  did  not  bifurcate 
in  the  next  section. 
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4.  Discussion  and  Conclusions 

Before  proceeding  with  further  discussion  of  our  numerical  results,  let 
us  try  to  lend  support  to  the  idea  that  these  tilted  and  C-shaped  plumes 
are,  in  fact,  seen  in  equatorial  spread  F.  We  present  experimental  radar 
backscatter  maps  of  meter  scale  plasma  irregularities  taken  during  equatorial 
spread  F.  Figure  11  shows  a  map  of  3  meter  irregularities  (provided  by  J.P. 
McClure)  using  the  Jicamarca  radar.  Similar  plots  can  be  found  in  Woodman 
and  La  Hoz  fl976]  .  Figure  12  shows  a  map  of  1  meter  irregularities  taken 
by  Tsunoda  [l98l]  using  the  ALTAIR  radar.  We  refer  the  reader  to  the  re¬ 
spective  papers  for  a  detailed  explanation  of  these  plots,  but  we  point  out 
that  the  Jicamarca  radar  scans  a  fixed  line  in  space,  and  plots  backscatter 
strength  as  a  function  of  time.  Therefore,  structures  caught  up  in  our 
postulated  plasma  shear  would  have  their  C-shaped  appearance  exaggerated  in 
the  Jicamarca  plots.  The  ALTAIR  radar,  however,  is  steerable;  and  it’s 
backscatter  plots  are  a  good  approximation  to  a  "snapshot"  of  the  backscatter 
strength  at  a  single  time.  In  both  plots  the  evidence  of  oppositely  tilting 
structures  at  high  and  low  altitudes  is  apparent.  In  making  comparisons  with 
these  small  scale  (^3m)  irregularity  radar  backscatter  maps  we  are  assuming 
that  these  maps  are  signatures  (Tsunoda,  1980)  of  the  large  scale  size  bubbles 
depicted  in  Figure  8.  That  is  the  steep  plasma  density  gradients  associated 
with  the  bubbles  in  Figure  8  drive  the  radar  backscatter  observed  irregulari¬ 
ties  .  The  westward  and  upward  motion  of  the  bubbles  depicted  in  Figure  8 
is  in  agreement  with  the  satellite  in  situ  measurements  of  McClure  et  al. , 

[1977]. 

It  is  our  belief  that  Che  arguments  advanced  in  this  paper  offer  the 
most  plausible  explanation  yet  of  the  qualitative  behavior  of  equatorial 
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Fig.  12  —  Maps  of  1  meter  irregularities  taken  by  the  steerable  ALT  AIR  radar  during 
equatorial  spread  F,  from  Tsunoda  [19811 .  These  maps  are  very  close  to  being  snap¬ 
shots  of  the  locus  of  irregularities  at  a  given  time.  Note  the  oppositely  tilting  struc¬ 
tures  at  different  altitudes  in  the  left-hand  map,  and  the  clear  bending  of  a  plume  into 
a  “C”  shape  about  an  altitude  of  about  400  km  in  the  right-hand  plot. 


spread  F  plumes:  "C"-shaped  structures  and  westward  plume  tilts  are  simply 
the  result  of  vertically  rising  spread  F  plumes  being  caught  up  In  the  ambient 
plasma  shear  both  as  they  rise  and  subsequently,  this  shear  being  the  natural 
consequence  of  a  neutral  wind  at  the  equator  and  E  regions  of  finite  conduc¬ 
tivity  connected  to  the  equatorial  F  region  along  magnetic  field  lines.  If 
there  Is  a  neutral  wind,  but  no  E  region,  the  plasma  will  move  at  the  wind  speed, 
ESF  bubbles  will  rise  vertically,  and  the  attendant  radar  hackscatter  maps 
will  show  non-tllted  (l.e. ,  vertical)  plumes,  as  exhibited  In  the  measurements 
of  Kelley  et  al.,  (1981).  Furthermore,  even  without  the  equatorial  F  region 
neutral  wind  the  numerical  simulations  show  that  E  region  Pedersen  conductiv¬ 
ity  can  have  a  dramatic  effect  on  ESF.  For  example,  the  results  of  section 
3  show  that  ESF  has  been  slowed  down  as  has  the  attendant  bubble  evolution. 

In  addition,  the  ESF  bubble  in  the  presence  of  an  E  region  Is  less  depleted 
than  without  the  E  region.  This  is  due  to  the  fact  that  the  E  region  has  a 
dramatic  effect  on  the  Induced  polarization  electric  field  (see  Fig.  9)  which 
causes  the  rise  of  the  bubble,  the  fringe  field  component  of  which 
(Zalesak  and  Ossakow,  1980)  determines  the  region  below  the  F  peak  from 
which  plasma  is  drawn  (i.e.,  which  plasma  makes  up  the  bubble).  The  pre¬ 
sence  or  absence  of  an  E  region  could  also  explain  why  bubbles  (with  large 
depletions)  stop  rising  at  altitudes  of  400-500  km  (see  McClure  et  al.,  1977) 
on  some  occasions,  but  not  on  others.  Indeed,  in  addition  to  the  height  of 
the  F  peak  and  bottomside  background  electron  density  gradient  scale  lengths 
(Ossakow  et  al.,  1979)  influencing  ESF  evolution,  the  E  region  conductivity 
could  also  determine  why  one  does  or  does  not  observe  ESF  even  when  the 
previously  mentioned  conditions  are  satisfied.  Thus,  the  E  region  (even  at 
night)  could  be  a  controlling  factor  in  the  formation  of  ESF  irregularities 


and  could  account  for  such  things  as  the  longitudinal  Influence  QBasu  and 
Kelley.  1977;  Livingston.  1980)  on  ESF  formation  and  phenomena. 

There  are  two  matters,  however,  which  bear  further  discussion:  1)  the 
question  of  why  the  plume  in  calculation  2LE  bifurcated  while  those  in  cal¬ 
culations  2L  and  2LEW  did  not;  and  (2)  the  question  of  where  along  the  edge  of 
the  primary  plumes  in  calculation  2LEW  one  should  expect  to  see  secondary 
instabilities.  If  we  note  that  an  equatorial  spread  F  plume  (bubble)  is 
nothing  more  than  an  inverse  plasma  cloud  "finger"  (i.e.,  an  elongated  region 
of  low  density  plasma  (an  ESF  "plume")  penetrating  a  region  of  high  density 
plasma  is  the  inverse  of  an  elongated  region  of  high  density  plasma  (a  plasma 
cloud  "finger")  penetrating  a  region  of  low  density  plasma),  we  find  that 
the  question  of  why  the  2LE  plume  bifurcated  and  the  2L  plume  did  not  has 
already  been  answered  for  us.  McDonald  et  al..  [l98l]  have  shown  in  their 
study  of  bifurcation  tendencies  of  plasma  cloud  fingers  that  the  critical 
quantity  determining  the  speed  with  which  a  plasma  finger  will  bifurcate  is 
M,  the  ratio  of  the  Pedersen  conductivities  inside  and  outside  the  finger. 
When  M  is  moderate,  in  the  range  2  to  10,  the  bifurcation  tendency  is  high, 
while  for  M  near  1  or  M  greater  than  100,  the  bifurcation  tendency  is  ex¬ 
tremely  small.  Looking  at  the  2LE  plume  (Fig.  7),  and  recalling  that  we 
have  a  background  Pedersen  conductivity  of  0.12  times  the  maximum  equatorial 
plane  Pedersen  conductivity  we  find  that  M~^  (the  relevant  quantity  since  we 
are  dealing  with  inverse  plasma  clouds  here)  in  the  2LE  plume  is  about  9, 
making  it  a  prime  candidate  for  bifurcation;  while  for  the  2L  plume 
(Fig.  6)  is  about  10‘*,  indicating  a  bifurcation  tendency  near  zero.  The 
question  of  why  the  2LEW  plume  did  not  bifurcate  is  a  little  harder  to 
answer.  Based  on  the  arguments  advanced  above,  the  2LEW  plume  would  have 
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been  just  as  likely  a  candidate  for  bifurcation  as  the  2LE  plume.  However, 
we  note  that  the  2LEW  plume  is  rising  into  a  region  of  very  strong  plasma 
shear.  Perkins  and  Doles  [l975l  have  shown  that  such  a  shear  would  provide 
a  stabilizing  mechanism  for  any  secondary  instability  (i.e.,  bifurcation) 
which  attempted  to  grow  on  the  topside  of  the  2LEW  plume,  although  the  geome¬ 
try  used  in  their  study  was  considerably  simpler  than  that  associated  with 
a  rising  ESF  plume.  We  advance  this  mechanism  as  a  plausible,  but  less  than 
totally  convincing, explanation  of  why  the  2LEW  plume  did  not  bifurcate,  only 
because  we  can  provide  no  other  at  this  time. 

The  question  of  secondary  instabilities  on  the  perimeters  or  in  the 
Interiors  of  the  plumes  is  in  many  ways  the  most  interesting  aspect  of  this 
study.  We  will  confine  ourselves  to  the  2LEW  plume  (Fig.  8)  at  late  time 
(2331  sec)  since  it  is  both  the  most  Interesting  and  the  most  realistic. 

Given  the  limited  spatial  resolution  of  these  numerical  studies,  it  must  be 
realized  that  the  actual  numerical  simulation  of  the  evolution  of  small-scale 
secondary  instabilities,  within  the  context  of  the  present  simulations,  is 
an  impossibility.  However,  we  ^  have  the  resolution  to  be  able  to  observe 
the  precursor  of  the  plasma  fluid  instabilities  that  we  believe  are  active; 
the  steepening  of  electron  density  gradients.  By  observing  the  location 
of  these  regions  of  steepening,  and  by  augmenting  this  procedure  with  a  cell 
by  cell  local  stability  analysis,  we  should  be  able  to  predict  both  the 
location  of  secondary  instabilities  and  the  mechanism  causing  them.  We  use 
the  term  "local  stability  analysis"  to  mean  a  local  evaluation,  numerical  in 
this  case,  of  the  generalized  gradient  drift  growth  rate  y  given  by 
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Note  that  ^  in  Eq.  (45)  above  includes  the  self-consistant  polarization 
electric  field  given  by  the  solution  to  the  potential  equation  (42).  The 
influence  of  this  term  on  y  is  large,  and  any  stability  analysis  which  were 
to  ignore  these  polarization  fields  would  rest  on  shaky  ground.  Among 
the  assumptions  implicit  in  the  application  of  Eq.  (45)  are:  1)  the  growth 
rates  y  are  large  compared  to  the  speed  with  which  the  primary  Rayleigh- 
Taylor  mode  is  evolving;  and  2)  the  k-vector  associated  with  the  growing 
perturbation  is  perpendicular  to  (which  gives  maximum  growth). 

Looking  at  Figure  8  at  the  latest  time  (2331  sec)  we  see  that  the  pri¬ 
mary  regions  of  steepening  are  two:  the  west  wall  of  the  plume  at  low  al¬ 
titudes  (below  370  km  for  this  particular  model)  and  the  east  wall  of  the 
plume  at  higher  altitudes  (above  370  km).  Local  stability  analysis  verifies 
that  these  are  precisely  the  regions  of  largest  growth  rate  for  the  gradient- 
drift/gravltational  Rayleigh-Taylor  Instability.  A  less  complete  analysis  of 
just  the  effect  of  an  eastward  neutral  wind  on  a  more  or  less  vertical  plume 
would  predict  that  only  the  west  wall  of  the  plume  would  be  unstable,  but 
this  analysis  neglects  the  effects  of  the  bending  of  the  plume  which  orients 
the  normally  stable  east  wall  of  the  plume  so  that  it  is  once  again  unstable 
to  the  gravitational  instability,  and  the  effects  of  the  polarization  electric 
field  produced  self-consistently  by  the  ionosphere-plume  system,  whose  effect 
is  to  mitigate  the  expected  wind-driven  Instability  over  most  of  the  plume 
structure.  We  conclude,  then,  that  for  this  particular  plume,  we  would  ex¬ 
pect  secondary  instabilities  along  the  west  wall  at  low  altitudes  and  along 
the  east  wall  at  high  altitudes,  with  the  "switch"  taking  place  at  about 
370  km  altitude.  If  these  Instabilities  eventually  cascade  down  to  smaller 
and  smaller  scale  sizes  (or  provide  the  steep  plasma  density  gradients 
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necessary  for  further  Instability  mechanisms),  eventually  reaching  the  1-3 
meter  scale  sizes  seen  on  backscatter  radar  maps,  we  would  expect  the  radar 
maps  to  trace  out  the  locus  of  the  west  wall  of  the  plume  at  low  altitudes, 
and  the  east  wall  at  high  altitudes,  giving  rise  to  an  even  more  exaggerated 
"C"  shape  than  that  of  the  simulation  plume  (bubble)  Itself  (see  the  ex¬ 
aggerated  "C"  traced  out  by  the  locus  of  steepened  gradients  in  Figure  8 
at  2331  sec). 

We  wish  to  close  this  section  by  briefly  reviewing  work  by  other  re¬ 
searchers  which  we  believe  has  relevance  to  the  results  presented  here. 

Two  recent  papers  have  shown  experimental  evidence  of  a  shear  in  east-west 
plasma  motion  in  the  equatorial  ionosphere:  Kudeki  et  al . .  [  19813  and 
Tsunoda  et  al. ,  [l98l]  .  Both  papers  show  evidence  of  an  Increase  in  eastward 
plasma  velocity  with  altitude,  in  agreement  with  the  behavior  we  postulate 
here  for  altitudes  below  the  peak  in  equatorial  F  region  Integrated  Pedersen 
conductivity.  It  is  our  belief  that  experimental  observations  at  even  higher 
altitudes  than  that  examined  in  the  above  papers  would  show  a  decrease  in 
eastward  plasma  velocity  with  altitude  at  these  higher  altitudes  (and 
possibly  even  westward  velocities),  in  a  manner  similar  to  that  shown  in 
Figure  3.  Both  Kudeki  et  al.,  [l98l3  and  Tsunoda  et  al.,Cl98l3  show 
evidence  of  a  plasma  velocity  reversal  point,  that  is,  an  altitude  below  the 
F2  peak  below  which  the  plasma  velocity  actually  becomes  westward  (as  the 
eastward  velocity  passes  through  zero).  The  simple  model  we  have  presented 
here  offers  no  explanation  for  this  phenomenon.  The  reason  is  that  we  have 
assumed  here  that  the  E  regions  connected  to  the  equatorial  F  region  along 
field  lines  are  passive  and  free  of  any  dynamics  of  their  own.  Actually, 
however,  these  E  regions  are  subject  to  strong  diurnal  tidal  neutral  winds. 
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which  are  westward  at  the  times  associated  with  spread  F.  At  equatorial  al¬ 
titudes  well  below  the  F2  peak,  the  field-line  Integrated  Pedersen  conductiv¬ 
ity  Is  dominated  by  that  of  the  E  regions,  and  hence  the  westward  neutral 
winds  In  the  E  regions  are  able  to  set  up  polarization  electric  fields 
which  Impress  themselves  on  the  equatorial  region,  causing  a  corresponding 
westward  drift  of  plasma  In  the  equatorial  plana  at  low  altitudes. 

In  fact,  had  we  retained  the  neutral  wind  terms  In  layers  1  and  3,  Eq. 
(38)  would  have  become 
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(46) 


where  U  ,  U  ,  and  U  are  the  east-west  neutral  winds  In  layers  1,  2, 

ni  n2  n3 

and  3  respectively.  If  U  Is  eastward  and  both  U  and  U  are  westward, 

n2  ni  ns  * 

It  Is  obvious  that  westward  plasma  velocities  will  exist  at  any  altitude 
for  v/hich 

Z  U  Z  U  , 

■  .E1.IL3 _ ElJi  (4 

U  '• 

n2 

This  effect  and  the  consequent  plasma  velocity  reversal  point  were  first 
described  by  Hellis  et  al. , [l 974],  whose  detailed  self-consistent  numerical 
model  of  the  E-and  F-reglon  neutral  gas  and  plasma  system  also  shows  both 
the  plasma  shear  and  an  altitude  at  which  the  eastward  plasma  velocity 
maximizes,  as  we  have  proposed  here.  In  fact,  at  extreme  equatorial  F 
region  altitudes,  plasma  well  away  from  the  equatorial  plane  (both  E  and  F 
region  plasma)  may  again  dominate  the  Integrated  Pedersen  conductivity,  and 
If  the  corresponding  neutral  winds  are  westward,  we  should  expect  to  see 


Z  < 
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westward  plasma  drifts  at  these  high  altitudes,  as  mentioned  above. 

The  Influence  of  contributions  to  field  line  Integrated  conductivity 
from  plasma  away  from  the  equatorial  plane  on  the  rise  of  ESF  bubbles  has 
also  been  addressed  by  Anderson  and  Haerendel  [l979].  In  their  model  they 
Incorporated  flux  tube  Integrated  quantities  of  electron  content  and  Pedersen 
conductivity  and  utilized  a  one- dimensional  sheet  model  for  the  bubble. 

This  latter  assumption  resulted  In  a  simple  algebraic  expression  for  the 
Induced  polarization  electric  field  inside  the  bubble  in  terms  of  the  flux 
tube  Integrated  quantities.  ESF  bubble  rise  in  the  collision  dominated 
Rayleigh- Taylor  regime  with  and  without  an  ambient  eastward  electric  field, 
was  investigated.  These  authors  noted  that  flux  tube  integrated  quanti¬ 
ties  could  have  a  considerable  influence  on  the  outcome  of  ESF  bubble  rise,  an 
observation  consistent  with  our  comparison  of  the  2L  and  2LE  cases.  Burke 
[19793  analyzed  the  effect  of  the  sunrise  "turning-on"  of  the  E  region  and  its 
subsequent  contribution  to  the  demise  of  ESF  bubbles  in  the  topside  F  region 
ionosphere  near  the  dawn  terminator.  A  simple  analytic  model  was  used  which 
showed  how  electric  fields  within  the  bubbles  could  be  discharged  through 
the  conducting  sunrise  E  region.  The  results  showed  that  the  conducting 
E  region  could  effectively  halt  the  upward  bubble  rise  velocity.  This  im¬ 
portance  of  an  E  region  is  consistent  with  our  findings  in  comparing  the  2L 


and  2LE  cases. 


5.  Future  Work 


The  "three  layer  model"  used  in  this  study  Is  simple  to  be  sure,  al¬ 
though  we  believe  It  adequately  describes  the  qualitative  behaviour  of  ESF 
plumes  In  terms  of  C-shaped  structures,  westward  tilts,  and  the  effects  of 
a  background  E  region.  Work  toward  Improving  the  model  and  Its  Input  Is 
ongoing  on  several  fronts.  First,  we  would  like  to  make  the  E  regions 
"active",  l.e.,  to  allow  external  forces,  such  as  neutral  winds  to  act  on 
the  E  region  plasma,  and  to  self-consistently  solve  the  continuity  equation 
there.  Second,  we  would  like  to  better  resolve  the  plasma  distribution  along 
magnetic  field  lines  by  adding  more  layers  to  the  model  to  represent  plasma 
between  the  equatorial  plane  and  the  E  regions.  The  total  number  of  layers 
might  be  seven  or  nine.  Third,  we  would  like  to  Incorporate  more  realistic 
models  of  electron  and  neutral  density  distributions,  external  electric  fields, 
neutral  winds,  and  chemistry  Into  the  models.  Sources  for  this  Information 
might  be  empirical  data  or  models  of  the  type  developed  by  Anderson  {l973] 
and  by  Forbes  and  Garrett  [19783.  Last,  but  certainly  not  least,  a  continuing 
effort  Is  being  made  to  keep  the  numerical  techniques  used  in  the  code  as 
close  to  state-of-the-art  as  possible.  A  recent  advance  [zalesak,  198 1] 
should  significantly  Improve  our  already  quite  good,  but  certainly  not 
perfect,  numerical  algorithms  for  solving  the  continuity  equation  In  the 


very  near  future. 
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01a  Attn  tooE  4700  T.  p.  Coffey  (25  as  if  un,  1  a  if  class) 
Ola  attn  Code  4701  JAa  D.  Brown 

01a  ATTN  Code  4780  Branch  head  (150  as  if  um,  1  cy  if  class) 

Ola  Attn  Code  7500 

01a  attn  Code  7550 

OlcY  attn  Code  7580 

01a  attn  Code  7551 

01a  ATTN  Code  7555 

01a  attn  Code  4730  t.  hcLean 

01a  attn  Code  4187 

COHMUOCft 

Navm.  Sea  Systems  Comham) 

WASHINGTCN,  D.C.  20362 

01a  attn  Capt  R.  Pitkin 

COMWOER 

naval  Space  Surveillance  System 
UAtl.GREN,  Va  22448 

01a  attn  Capt  j.  h.  Burton 

Cfficer-In-Chahge 
Navm.  Surface  weapons  Center 
VAiiTE  Oak.  Silver  Spring,  it)  20910 
01a  Attn  Code  F31 

Director 

Strategic  Systems  Project  Office 
Department  of  d*  Navy 
Washington,  D.C.  20376 
01a  ATTN  NSP-2141 
01a  ATTN  NSSP-2722  Fred  Wimberly 

ComucER 

Naval  Surface  weafons  Center 
Dah-Oien  laboratory 
DAtCGREN.  Va  22448 

OViY  ATTN  Code  Uf-14  r.  Butler 

OFFICE  OF  naval  RESEARCH 
Arlington,  va  22217 
01a  ATTN  Code  465 
Ola  attn  Code  461 
Ola  attn  Code  402 
Ola  attn  Code  420 
Ola  attn  Code  421 

COMWOER 

aerospaq  defense  Comwwo/DC 
Department  of  the  air  Force 
ENT  AFB.  CO  80912 

01a  attn  DC  t*.  LONG 

COWWCCfl 

Aerospace  Defense  Commuo/xpO 

0EPART»€NT  OF  THE  AlR  FORCE 

ENT  AFB.  CO  80912 
01a  ATTN  XPDOO 

01a  Attn  xp 

Air  Force  Oeophysics  laboratory 
hanscom  AFB,  HA  01731 

01a  ATTN  OPR  HAROLD  GARDNER 

01a  attn  opr-1  jaws  C.  Xwick 

01a  attn  ucb  kewcth  s.  w.  Chmvion 
Ola  attn  opr  /4.VA  T.  Stair 
01a  attn  PHP  Jules  Aarons 
01a  Attn  PHD  Jurgen  Buchmj 
01a  attn  PHD  John  p.  hullen 
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If  VCN>ONS  UBORATORY 

KIRTVANO  AFB.  NH  87117 
Ola  ATTN  SUL 

Ola  ATTN  CA  ArTH*  H.  GUENIWR 
OICY  ATTN  NTYC  ILT  G.  KRAJCl 


^tac 

Patrioc  AFB.  Pl  829E5 
Ola  AHN  TF/11AJ  WIlfT 
01a  ATTN  TH 

AIR  Force  avionics  LAaoRATORY_ 

'i«IGHT-PATrERSON  AfB.  OH  XSRJJ 

Ola  Attn  AAC  haee  hunt 
Ola  Attn  AAO  Allen  joiwson 

Depotv  Chief  of  Staff 

KESEARCM,  DEVELOFNENT,  i  ACQ 
DERNTTMENT  OF  THE  AlR  FORCE 
WASHINCTON.  D.C.  20550 
Ola  Attn  AFRDQ 

HEAOaUARTERS 

aECTRONIC  SVSTERS  DIVISION/XR 

Defwtment  of  the  Air  Force 
Hanscon  AFB-  HA  01731 
Ola  Attn  XR  J.  L£AS 

HEADQUARTERS  „ 

cLECTROHic  Systems  Oivisiqn/YSEA 
Defnitment  of  t«  Air  FoRa 
HANSCON  AFB.  HA 01732 
Ola  ATTN  YSEA 

HEADQUARTERS 

Electronic  Systems  jivision/uc 
[jep/wtment  of  the  air  Force 
HANSCON  AFB,  HA  01731 

Ola  ATTN  WX  maj  j.l.  Clark 

CQNNAICER 

Foreion  Teoiwlosy  i/IVISIOW,  aFSC 
miONT-PATTERSON  AFB,  UH  MSA53 
Ola  ATTN  NlCD  LIBRARY 
Ola  ATTN  ttOH  B.  BALLARD 

NONE  AIR  UEVELOFMENT  CENTER.  aFSC 
tRIFFISS  AFB.  NY  13A41 

Ola  ATTN  UOC  LIBRARY/TSLU 
Ola  Attn  ocSE  v.  Coyne 

S*SO/SZ 

Post  office  box  92960 
MDRLWAY  Postal  Center 
Los  Angeles.  CA  90009 
(SPAa  Defense  Systems) 

Ola  Attn  SZJ 

STRAnGIC  air  CONMAAO/XPFS 
Offutt  AFB.  w  08II3 

Ola  ATTN  XPFS  Naj  B.  Stephan 
Ola  ATTN  ADHATE  Haj  Bruce  Bauer 
Ola  ATTN  NfiT 

Ola  ATTN  DOK  Chief  Scientist 
SAMSO/SK 

?.  0.  Box  92960  , 

HORLONAY  Postal  Cehtb 
LOS  Angles.  Ca  90006 

Ola  Attn  SKA  (Sp/ce  Coat  Systems)  m.  Clavin 


Department  of  tNEna 
Albuquerque  Operations  Office 
P.  0.  Box  SROO 
ICBUQUERQUE.  141  87115 

Ola  Attn  doc  con  for  D.  Shernooo 

Department  of  tNERa 
LIBRARY  Room  G-0>(2 
mashington.  D.C.  205R5 

Ola  attn  doc  con  for  a.  labonitz 

£686,  Inc. 

LOS  Alamos  Division 
P.  0.  BOX  809 
LOS  /R-AMOS.  NH  855RR 

Ola  attn  doc  con  for  j.  Bieedlove 

UNIVERSITY  OF  CALIFORNIA 
LAMRENCE  LIVERMORE  LABORATORY 
P.  0.  BOX  808 
LIVERMORE,  CA  9M550 

Ola  attn  doc  con  for  Tech  Info  Deft 
Ola  attn  DOC  CON  FOR  L-389  R.  On 
Ola  attn  DX  CON  FOR  L-31  R.  HAGER 
Ola  attn  doc  con  for  L-R6  F.  Semard 

LOS  Alamos  Scientific  laboratory 

P.  0.  BOX  1663 

Los  Al>*<js.  1*1  875*6 

Ola  attn  Doc  Con  for  j.  holcott 
Ola  attn  Doc  Con  for  R.  f.  Taschbi 
Ola  attn  Doc  Con  for  t.  x«€s 
Ola  ATTN  Doc  Con  for  J.  Halix 
Ola  attn  Doc  Con  FOR  K.  jeffries 
Ola  MTTN  Doc  Con  FOR  J.  2iw 
Ola  attn  Doc  Con  for  P.  Keaton 
Ola  attn  doc  Con  for  D.  i«stervelt 

Sam)ia  Laboratories 
P.  0.  BOX  5800 
Albuojerque.  nh  87115 

01a  attn  Doc  Con  for  j.  Martin 
01a  Attn  Doc  Con  for  w.  i*o»< 

01a  Attn  Doc  Con  for  a.  Thornbrough 
01a  Attn  Doc  Oon  for  f.  t«i«r 
01a  Attn  Doc  Con  for  D.  Dah.gr£n 
01a  ATTN  Doc  Con  for  31R1 
Ola  ATTN  Coc  Con  for  space  PRooEa  uiv 

Sandia  laboratories 

LIVERMORE  LABORATORY 
P.  0.  BOX  969 
LIVERMORE,  CA  9<650 

Ola  Attn  Doc  Con  for  b.  hjii»»€y 
Ola  attn  doc  lon  for  t.  Co« 

OFFICE  OF  Military  Application 
Department  of  tNEna  ^ 

WASHINGTON,  D.C.  206*6 

Ola  ATTN  Loc  Con  for  D.  bale 

OTHER  OOVERNttNT 

Centrn.  Intelligence  AcENa 
ATTN  KQ/Sl,  KM  5G**8.  HO  BLDG 
WASHINGTON.  D.C.  20605 

01a  Attn  OSI/PSID  Rm  5F  19 

Department  of  ConcRCE 
national  Bureau  of  Staniards 
washingtcn.  D.C.  2023** 

(All  Corres:  attn  Sec  Officer  for) 

01a  attn  R.  MOORE 


hS^nAFB.  Ca  92*09 
(HIHIT^) 

Ola  Attn  MHH. 


lTC  KEfWEOY 


COMMfDCR 

RofC  air  Deveuof*nt  Center,  AfSC 
HAABCOM  AFB,  MA  01731 

Ola  attncEP  a.  lorentzen 


Institute  for  Telecom  Sciemxs 
nationm.  Telecommunications  8  Info  Admin 
bouloer.  Co  80303 

Ola  ATTN  A.  JEAN  (UNCLASS  ONLY) 
Ola  ATTN  H.  OTLAUr 
Ola  ATTN  u.  Crohbie 
Ola  ATTN  L.  BERRY 
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national  oceanic  i  atmospheric  ^^in 
tNVIROtWENTAL  KESEARCH  LABORATORIES 
APARTMENT  OF  COM«RC£ 
fiOOLIXR.  Co  110302 

Ola  ATTN  R.  Grubb 

Ola  ATTN  AERONOMY  LAB  G.  KEID 

iJEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  Corporation 
P.  0.  Box  92957 
Los  Anseus,  Ca  90009 
Ola  Attn  I.  Garfunkel 
Ola  Attn  T.  Salmi 
Ola  Attn  v.  josepmson 
Ola  Attn  S,  Boner 
Ola  ATTN  N.  STOOCNELL 
OfcY  Attn  D.  Olsen 


ANJ^YTICAL  SySTWS  ENGIKERINS  CORP 

s  Old  Concord  koao 
BURLINGTON.  IV  01803 

Ola  ATTN  Radio  Sciences 

Associates.  Inc. 

BERkLev.  ca  9<t701 

01a  Attn  j.  workman 

Boeing  CDm>any,  The 
P.  0.  Box  3707 
Seattle,  ha  9812A 

Ola  ATTN  G.  Keister 
gia  Attn  u.  purray 
Ola  ATTN  6.  HAa 
Ola  Attn  j.  kemcy 


Gewral  Electric  Cowany 
SPACE  Division 
VALLEY  Forge  Space  Center 
God^d  a.TO  King  of  Prussia 
P.  0.  Box  8555 
Philadelphia,  Pa  19101 

Ola  Attn  h.  h.  Bortner  Sp/ce  Sci  lab 

Geicral  Electric  Co»*>any 
P.  0.  Box  1122 
Syracuse,  nY  13201 

Cla  mttn  f.  reibert 

^NERAL  aECTRIC  COfPANY 
lEffO-CENTER  FOR  ADVANCED  STUDIES 

816  State  Street  (P.O.  dramer  m 
Santa  Barbara,  Ca  93102 
Ola  Attn  iiasiac 
gia  Attn  Don  Chaicur 
Ola  Attn  Tom  Barrett 
Ula  ATTN  Tim  Stephans 
Ola  ATTN  HARREN  S.  AN/MN' 

Ola  ATTN  WILLIAM  MCNAMMA 
Ula  Attn  b.  ixucia 
Ola  ATTN  MAtt  STANTON 

general  aECTRIC  TECH  SERVICES  CO.,  iNC. 

Court  Street 
SYRACUSE,  NY  13201 

Ola  ATTN  6.  HILLMAN 


California  at  San  Diego,  lniv  of 

P.O.  aox6(M9 

San  Diego.  CA  92106 


iXMl  ENGIfEERINS  COWANY,  InC. 
CitwiNGS  Research  Park 
HUNTSVILLE,  AL  35807 

Ola  ATTN  Romeo  a.  Deliberis 


Laboratory,  nc. 

555  TEotioLOGY  Square 
Camkiose,  Na  02139 
Ola  Attn  D.  B,  Cox 
Ola  Attn  j.  P.  Gilmore 


general  research  Corporation 
Santa  Barbara  Division 
P.  0.  BOX  6770 
Santa  Barbara.  Ca  93111 
gia  Attn  joHt  ise  Jr 
01a  Attn  joel  Garbwino 


University  of  Alasu 
Fairbaws,  ak  99701 

A,  Seohity  officer) 

gia  Attn  t.  h.  Davis  (Uha  only) 

OICY  Attn  Neal  Brom  (UNa  Only) 

Ola  Attn  Teowicin.  library 


COnSAT  Laboratories 
lintnicun  Road 
Clarksbirg,  It)  2073M 
Ola  Attn  G.  myoe 

CORieu  UNIVERSITY 

Ola  MTTN  0.  r.  Farley  Jr 

CLECT^ACE  SySTENS,  InC. 
oox  1359 

Richardson.  IK  75080 
Ola  Attn  h.  logston 
Ola  ATTN  Security  (Paul  Phillips) 


£SL  Inc. 


m96  Java  Drive 


SownvALE,  Ca  3R086 
Ola  Attn  j.  Roberts 

Ola  ATTN  J/*€S  IVRSM«X 
Ola  ATTN  C.  W.  PRETTIE 


G1I  Sylvania,  Inc. 

^ECTRON^  Systems  G«P-Eastern  Div 

neewah,  ha  Q219A 

Ola  Attn  marshal  Cross 


lULinuia,  univtRSJTY  OF 

OOBU  HAU 
ISO  DAVEIf>ORT  HOUSE 
CHAAPAIGN,  IL  61820 

gj.L  CORRES  Attn  Dan  rtaELLAic) 
Ola  For  k.  Yen 


IjernuTE  for  defense  analyses 
MOO  army-navy  Biive 


Arlington,  Va  22202 
Ola  Attn  J,  n,  aein 
gia  Attn  Ernest  bauer 
Ola  Attn  HANS  hdlfhard 

Ula  ATTN  JOEL  BENGSTDN 
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HSS.  INC. 

I  Alfred  Circu_ 
liECFQRO.  Na  01730 

01CY  Attn  Domld  hansen 

Iwn.  Tel  &  Telesrafh  Corporation 
SOO  Washinston  Anenue 
NUTLEY,  «J  07110 

OICY  Attn  Technical  lierary 

JVYCOR 

1M01  Camno  Del  Nar 
Del  Nar.  Ca  920W  ^ 

Ola  Attn  S.  R.  Golomn 

JOHNS  Hopkins  University 
Applied  Physics  Laboratory 
joms  Hopkins  Road 
LAUREL.  Ho  20B10 

Ola  Attn  Document  Librarian 
Ola  Attn  Thomas  Pote*a 
Ola  Attn  Joiti  Dassoulas 

LocxtcED  Missiles  S  Space  Co  Inc 
P.  0.  Box  504  _ 

auNKYVAU.  Ca  94088 

Ola  Attn  Dept  60-12 
Ola  Attn  0.  R.  Churchiu 

LogocEO  Missiles  ami  Space  Co  Inc 
3251  HANOVER  street 
Palo  alto,  u  94304 

Ola  ATTN  Martin  Halt  Dept  52-10 
Ola  ATTN  KIOWRO  6.  J0M6QN  UEPT  52*12 
Ola  Attn  h.  l.  Iw  Dept  52-12 

aaman  sciences  Corp 
P.  0.  sox  7463 
COLORADO  Springs,  lo  a0!63 
Ola  ATTN  T.  HEAGHER 

LINUBIT  CORP 
10453  moselle 
SAN  Diego.  Ca  92121 

Ola  ATTN  IRMIN  Jacobs 

n.l.T.  LINCOLN  LABORATORY 

P.  0.  Box  73 
LEXIfBTON.  MA  02173 

Ola  ATTN  David  n.  To«a.e 

Ola  ATTN  P.  HALDRON 

Ola  Attn  l.  loug«.in 
Ola  ATTN  D.  Clark 

MARTIN  IHRIETTA  CoRP 
URLANDO  DIVISION 
P.  0.  BOX  5837  „ 

ORLMCO.  FL  3305 

Ola  Arm  k.  meffner 

PWSICM.  DYNAMICS  INC. 

P.  0.  Box  3027 
BELLEVUE.  HA  98009 

Ola  Aim  E.  J.  Fremoiw 

Physicn.  dynamics  Inc. 

P.  0.  box  1(S& 

OAKUM).  Ca.  94610 
ATTN:  A.  Thomson 

R  i  D  ASSOCl^S 
p.  0.  BOX  9695 
(Mina  Del  Rey.  Ca  90291 

01a  ATTN  Forrest  Gilmore 
Ola  Arm  uiyni  Gabbard 
Ola  Arm  william  B.  might  jr 
Ola  Arm  Robert  f.  ulevier 
Ola  Arm  hiuiam  j.  karzas 
Ola  Atm  H.  Ory 
Ola  Arm  C.  MacDonald 
Ola  Arm  R.  Turco 


RAN)  Corporation.  The 
1700  Main  Street 
Santa  Monica.  Ca  90406, 

Ola  Arm  Cullen  Crain 
Ola  Arm  Eo  Becrozian 

Riverside  Research  Ihstitute 
SO  HEST  End  aivenk 
NEH  York,  ny  10^ 

Ola  Arm  vince  Trapani 

SCIENQ  /miCATIONS.  iNC. 

P.  u.  BOX  2351  _ 

la  JOLLA.  Ca  92038 

Ola  Arm  uws  m.  linson 
Ola  Arm  DANia  a.  man.in 
Ola  ATTN  D.  SAOe 
Ola  Arm  t.  a.  straker 
Ola  Arm  Curtis  a.  SMim 
Ola  Arm  jacx  mcDougall 

raytnecn  Co. 

528  Boston  Post  Road 
SUDBURY,  Ha  01776 

Ola  Arm  Barbara  Adams 

SCIENCE  M>FL1UTI0NS.  iNC. 

1710  GOODtIDGE  DR. 

MdEAN,  VA  22102 

Arm:  J,  Cockayne 

Lockheed  Missile  &  Sp^  Co.,  Inc. 

HUNTSVILLE  RESEARCH  &  EhGR.  CTR. 

4800  Btacford  Drive 
Huntsville.  Alabama  35807 
Arm:  u*u  H.  Davis 

ncDoMca  DOUGLAS  Corporation 
5301  HOLSA  AVENUE  „  , 

HUNTINGTON  UEAOI.  CA  92647 
Ola  Arm  N.  Harris 
Ola  Arm  J.  hcxa-e 
Ola  ATm  UEORGE  moz 
Ola  ATm  N.  Olson 
Ola  AHN  R.  H.  iVUJHIIN 
Ola  Arm  Technical  Library  services 

MISSION  Research  Corporation 
735  State  Street 
Santa  mbara,  ca  93101 
Ola  Arm  P.  Fischer 
Ola  Arm  H.  F.  Crevier 
Ola  Aim  STEVEN  L.  UUTSCHE 
Ola  ATTN  D.  SimNIELD 
Ola  Arm  R.  dogusch 
Ola  Arm  R.  menricx 
Ola  Arm  Ralph  Kilb 
Ola  Arm  Dave  Sokle 
Ola  Arm  F.  Fajen 
Ola  Arm  H.  Scheibe 
Ola  Aim  CowAO  L.  longhire 
Ola  Arm  hwren  a.  Socueter 

Mitre  Corporation.  The 
P.  0.  Box  20B 
BEDPORD,  ha  01730 

Ola  Arm  jom  Hdrgawtfw 
Ola  Arm  6.  harddc 
Ola  Aim  C.  t.  Callahan 

MITRE  Corp 

HESTGATE  RESEARCH  PARK 

1820  Douy  Madison  Blvd 
ncUAN,  VA  22101 

Ola  Arm  h.  mau 
Ola  Arm  h.  Foster 

Pacific-sierra  Researo*  Corp 
1456  Cloverfielo  a.vD. 

SANTA  MONICA.  CA  90404 

Ola  Arm  £.  c.  Fiao  jr 
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PEMtsYLVMiA  state  University 

tONOSPMERE  RESEMCH  LAB 

318  aECTRICAL  ENGIICERINB  LAST 

UNIVERSITY  Park,  PA  16802 

(No  OASSIFIED  TO  THIS  AOCRESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LAB 

Photoactrics,  Inc. 

4<C  Mmett  Road 
lexinston,  Na  02173 

OICY  Attn  Irving  L.  Kofsky 

TEOttOLOGY  International  Corp 
75  WIGGINS  Avenue 
Beoforo.  IV  01730 

01a  Attn  w.  p.  Boouist 

TMW  Defense  t  Space  Svs  Sucuf 

One  Space  Park 

KEOONDO  Beach,  U  50278 

01a  ATTN  R.  K.  ft.£BUCM 
01a  Attn  S.  «ltsq«.eb 
01a  Attn  u.  Dee 


SRI  International 
333  kavenswood  Avenue 
nEN.0  Park,  U  9N02S 


01a 

Ola 


Ola 

01a 


Attn  Donald  neilson 
Attn  Alan  Burns 


01a  ATTN  6.  Snith 
Ola  Attn  l.  l.  Cobb 
Ola  Attn  David  A.  Johnson 
01a  Attn  Walter  G.  Chesnut 
Ola  attn  Charles  l.  Rino 
01a  Attn  Walter  jaye 
Ola  Attn  n.  Baron 
Ola  Attn  Ray  l.  leaoabraw 
Ola  Attn  6.  Carpenter 
01a  Attn  G.  Price 
01a  attn  j.  Peterson 


Attn  R.  hmce,  J«. 
attn  v.  Gonzales 


0la  attn  D.  HcDaniel 
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IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
UNCLASSIFIED  ONLV 


Pitast  district!  ona  copy  to  Mch  of  tho  follOHlnq  pooplo: 


Naval  Rtsoarcti  Laboratory 
Uasb1n9t00t  O.C.  2037S 

Or.  P.  Niopt  •  £odt  4101 

Dr.  N.  %lor  .  Codt  4141 

Or.  C.  SntacMwIci  -  Coda  4187 

Or.  J.  teodMO  •  Coda  4180 

Sclanca  Applications.  Inc. 

12S0  Prospact  Plaza 
La  Jolla.  California  02037 

Or.  0.  4.  Hmlitt 
Dr.  L.  Linton 
Or.  0.  Sachs 

Oiroctor  of  Spaca  and  Cnvlronoantal 
Laboratory.  NQAA 
Bouldar.  Colorado  00302 

Or.  A.  filann  Jaan 
Or.  0.  M.  Ada« 

Or.  0.  N.  Andarson 
Or.  K.  Oavlas 
Or.  A.  f.  Oonnally 

A.  f.  Gaaphysics  Laboratory 
L.  G.  Hansoa  Plaid 
Oadford,  Mass.  01730 

Or.  T.  Clkint 
Or.  M.  Swidar 
Mrs.  A.  Sapalyn 
Or.  J.  N.  Porbas 
Or.  r.  J.  Kanashaa 
Or.  J.  Aarons 

Offica  of  Naval  Aasaarch 
000  North  Quincy  Stroat 
Arllofton.  Virpinia  22217 

Or.  H.  Nullanay 

CdMOdtr 

Naval  Cloctronics  Laboratory  Cantor 
San  01090.  California  921S2 

Hr.  A.  Aoso  -  Coda  S321 

U.S.  Amy  Abardaan  Attaarch  and 
Oavalopaant  Cantor 
Ballistic  Nsatrch  Laboratory 
Abardaac.  *0 

Or.  J.  Halatrl 


CoMindar 

Naval  Air  SystoBS  Coanand 
O^rtaant  of  tN  Navy 
Vashlnyton,  O.C.  20360 

Or.  T.  Czuba 

Harvard  Univarsity 
Harvard  Souara 
CaHbrfdda,  Hass.  02130 

Or.  N.  B.  HcElroy 
Or.  A.  Lindzan 

Pannsylvanla  Stata  UnIvarsUy 
UnIvarsUy  Part,  Pannsylvanla  16802 

Or.  J.  S.  Nisbat 
Or.  P.  A.  AohrbavBh 
Or.  L.  A.  Carpantar 
Dr.  N.  Lao 
Or.  A.  Olvany 
Or.  P«  Bannatt 
Or.  C.  Klavans 


UnIvarsUy  of  CallfomU,  Barfcalay 
Bartalay,  California  9472Q 

Or.  N.  Hudson 

Utah  Stata  Univarsity 
ath  N.  and  6th  Stroots 
Lopan,  Utah  84322 

Or.  P.  N.  Banks 
Dr.  A.  Harris 
Dr.  K.  Bakar 

Comall  UnIvarsUy 
Ithaca.  Haw  fort  148S0 

Dr.  H.  E.  SiMrtz 
Dr.  A.  Sudan 
Or.  0.  Parlay 
Dr.  H.  Kallay 


NASA 

Goddard  Spaca  Pllpht  Cantor 
Graanbalt,  NarvUnd  2077] 

Or.  S.  Chandra 
Or.  X.  Naada 
Or.  A.P.  Sanson 

Princaton  UnIvarsUy 
Plasaa  Physics  Laboratory 
Princaton.  Na«  Jtrsay  08540 

Dr.  P.  Parkins 
Dr.  E.  Prlasan 

tnstUuta  for  Dafansa  Analysis 
400  Anay/Navy  Driva 
ArMnpton.  Vlrslnla  22202 

Or,  C.  Bauar 

UnIvarsUy  of  Maryland 
Collapa  Park.  «  20742 

Dr.  K.  Papadepeulos 
Dr.  £.  Ott 

UnIvarsUy  of  PUtsburoh 
Pittsburgh,  Pa.  1I>2d1 

Or.  N.  Zabusky 
Dr.  N.  Olondl 

Dafansa  Oocmiantatlon  Canttr 
Caaaron  Station 
ATaxandrfa.  M  22314 

(12  coplas  If  opan  publication 
othanvisa  2  coplas) 

12CY  Attn  TC 

UnIvarsUy  of  California 

Los  Alaaos  Sclantific  Laboratory 

J.IO.  HS.664 

Los  Alaaos,  Now  HaiUo  87545 

H.  Ponqract 
0.  Sfaons 
6.  larasch 
L.  Duncan 

Nassachusatts  InstUuta  of  Tachnelogy 
Plasai  Puslon  Cantor 
Library,  NH16-262 
Ciabrldgt.  HA  02139 


Univarsity  of  California,  Los  Angolas 
SOS  Hlllgard  Avanua 
os  Angafas.  California  90024 

Or.  P.  V.  Coronltl 
Or.  C.  Kannal 
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